% for beta = 0.00:0.05:0.60
%     for c_0 = 1:0.1:0.1*ceil(10*2*(1+beta)/(1-beta))
%         % clear;
%         clear_temp_folder
% 
%         N = 128;
% 
%         data = struct();
%         parfor j=1:N
%             [data(j).data, data(j).l1_norm, data(j).convergence, data(j).id, data(j).k_end] = runparallelcode(j,beta,c_0);
%         end
% 
%         clear N
% 
%         % save("./matfiles/"+ string(datestr(datetime(),'mmddHHMM')))
% 
% 
%         save("./matfiles/"+ sprintf('beta_%.2f_c0_%.1f.mat', beta, c_0));
% 
% 
%         % play_sound
% 
%     end
% end


for beta = 0.00:0.05:0.60
    % for gamma = 0.8:0.01:0.99
    for gamma = 0.991:0.001:0.999
        % clear;
        clear_temp_folder

        N = 1024;
 
        data = struct();
        parfor j=1:N
            [data(j).data, data(j).l1_norm, data(j).convergence, data(j).id, data(j).k_end] = runparallelcode(j,beta,gamma);
        end

        clear N

        % save("./matfiles/"+ string(datestr(datetime(),'mmddHHMM')))


        save("./matfiles/"+ sprintf('beta_%.2f_gamma_%.3f.mat', beta, gamma));


        % play_sound

    end
end




% for beta = 0.55:0.05:0.7
%     beta
% end
% 
% for gamma = 0.9:0.01:0.99
%     gamma
% end